Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis.

The multidrug and toxin efflux (MATE) family participates in numerous biological processes and plays important roles in abiotic stress responses. However, information about the MATE family genes in Torreya grandis remains unclear. In this study, our genome-wide investigation identified ninety MATE genes in Torreya grandis, which were divided into five evolutionary clades. TgMATE family members are located on eleven chromosomes, and a total of thirty TgMATEs exist in tandem duplication. The promoter analysis showed that most TgMATEs contain the cis-regulatory elements associated with stress and hormonal responses. In addition, we discovered that most TgMATE genes responded to abiotic stresses (aluminum, drought, high temperatures, and low temperatures). Weighted correlation network analysis showed that 147 candidate transcription factor genes regulated the expression of 14 TgMATE genes, and it was verified through a double-luciferase assay. Overall, our findings offer valuable information for the characterization of the TgMATE gene mechanism in responding to abiotic stress and exhibit promising prospects for the stress tolerance breeding of Torreya grandis.

Genome and Transcriptome Analysis of the Torreya grandis WRKY Gene Family during Seed Development.

Torreya grandis, an economically significant evergreen tree species exclusive to subtropical China, is highly valued for its seeds. However, the seed development process of T. grandis remains relatively unexplored. Given the pivotal role WRKY transcription factors (TFs) play in coordinating diverse cellular and biological activities, as well as crucial signaling pathways essential for plant growth and development, and the lack of comprehensive investigation into their specific functions in T. grandis, our study investigated its genome and successfully isolated 78 WRKY genes and categorized them into three distinct clades. A conserved motif analysis unveiled the presence of the characteristic WRKY domain in each identified TgWRKY protein. The examination of gene structures revealed variable numbers of introns (ranging from zero to eight) and exons (ranging from one to nine) among TgWRKY genes. A chromosomal distribution analysis demonstrated the presence of TgWRKY across eight chromosomes in T. grandis. Tissue-specific expression profiling unveiled distinctive patterns of these 78 TgWRKY genes across various tissues. Remarkably, a co-expression analysis integrating RNA-seq data and morphological assessments pinpointed the pronounced expression of TgWRKY25 during the developmental stages of T. grandis seeds. Moreover, a KEGG enrichment analysis, focusing on genes correlated with TgWRKY25 expression, suggested its potential involvement in processes such as protein processing in the endoplasmic reticulum, starch, and sucrose metabolism, thereby modulating seed development in T. grandis. These findings not only underscore the pivotal role of WRKY genes in T. grandis seed development but also pave the way for innovative breeding strategies.

Advancements in the Cultivation, Active Components, and Pharmacological Activities of Taxus mairei.

Taxus mairei (Lemée and H.Lév.) S.Y.Hu, indigenous to the southern regions of China, is an evergreen tree belonging to the genus Taxus of the Taxaceae family. Owing to its content of various bioactive compounds, it exhibits multiple pharmacological activities and has been widely applied in clinical medicine. This article comprehensively discusses the current state of cultivation, chemical constituents, applications in the pharmaceutical field, and the challenges faced by T. mairei. The paper begins by detailing the ecological distribution of T. mairei, aiming to provide an in-depth understanding of its origin and cultivation overview. In terms of chemical composition, the article thoroughly summarizes the extracts and monomeric components of T. mairei, unveiling their pharmacological activities and elucidating the mechanisms of action based on the latest scientific research, as well as their potential as lead compounds in new drug development. The article also addresses the challenges in the T. mairei research, such as the difficulties in extracting and synthesizing active components and the need for sustainable utilization strategies. In summary, T. mairei is a rare species important for biodiversity conservation and demonstrates significant research and application potential in drug development and disease treatment.

Diterpenoids from Torreya grandis and their cytotoxic activities.

Eight previously undescribed diterpenoids, along with eleven previously reported analogues, were obtained from the supercritical CO2 extracts of Torreya grandis aril. The structures of these compounds were elucidated based on HRESIMS, NMR, ECD, and single-crystal X-ray diffraction data. In the MTT assay, compound 18 exhibited significant inhibitory effects on two human colon cancer cell lines, HT-29 and HCT 116 cells, with IC50 values of 7.37 µM and 6.55 µM, respectively. It was found that compound 18 induced apoptosis and significantly inhibited the migration of HCT 116 colon cancer cells in a concentration-dependent manner.

TgLUT1 regulated by TgWRKY10 enhances the tolerance of Torreya grandis to drought stress.

Drought stress is a major abiotic stress which severely reduces the plant growth and limits agricultural productivity. Previous studies have demonstrated that lutein directly synthesized by the carotenoid epsilon-ring hydroxylase gene (LUT1) played crucial roles in regulating drought response. Notwithstanding the myriad studies on LUT1's response to drought stress in certain plant species such as Arabidopsis, the precise function mechanisms within tree species remain ambiguously understood. Our study reveals that under drought stress, TgLUT1, a novel LUT gene instrumental in ß-lutein biosynthesis, was markedly up-regulated in Torreya grandis. Subcellular localization assay indicated that TgLUT1 protein was localized to chloroplasts. Phenotypic analysis showed that overexpression of TgLUT1 enhanced the tolerance of tomato to drought stress. Overexpressing of TgLUT1 increased the values of maximal photochemical efficiency of photosystem II (Fv/Fm), net photosynthetic rate (Pn) and non-photochemical quenching (NPQ), and reduced the accumulation of hydrogen peroxide (H2O2), malondialdehyde (MDA) content and electrolyte leakage percentage in response to drought stress. Furthermore, overexpression of TgLUT1 decreased the stomatal conductance to reduce the water loss rate exposed to drought stress. In addition, yeast one-hybrid assay, dual luciferase assay system and qRT-PCR results showed that TgWRKY10 down-regulated by drought stress inhibited the expression of TgLUT1 by directly binding to the TgLUT1 promoter. Collectively, our results show that TgWRKY10, down-regulated by drought stress, negatively regulates the expression of TgLUT1 to modulate the drought stress response. This study contributes to a more comprehensive understanding of LUT1's function in the stress responses of economically significant forest plants.

Metabolomics analysis of the nutraceutical diversity and physiological quality of Torreya yunnanensis seeds during cold storage.

This study investigated how cold storage affects the nutraceutical diversity and physiological quality of Torreya yunnanensis seeds, using a widely targeted UPLC-MS/MS-based metabolomics analysis. The 373 identified metabolites were divided into nine categories: lipids, phenolic acids, amino acids and derivatives, organic acids, nucleotides, saccharides, vitamins and alcohols. Among them, 49 metabolites showed significant changes after 3 months of cold storage, affecting 28 metabolic pathways. The content of amino acid-related metabolites significantly increased, while the content of sugar-related metabolites decreased during storage. Notably, the content of proline acid, shikimic acid, α-linolenic acid and branched-chain amino acids showed significant changes, indicating their potential role in seed storage. This study deepens our understanding of the nutraceutical diversity and physiological quality of T. yunnanensis seeds during storage, providing insight for conservation efforts and habitat restoration.

Integrated Metabolomics, Transcriptome and Functional Analysis Reveal Key Genes Are Involved in Tree Age-Induced Amino Acid Accumulation in Torreya grandis Nuts.

Torreya grandis is native Chinese tree species of economic significance, renowned for its long lifespan and the rich nutritional value of its nuts. In this study, we analyzed the morphological characteristics, metabolites, associated gene expressions, and regulatory mechanism in nuts from young (10 years old) and old (1000 years old) T. grandis trees. We observed that the length, width, and weight of nuts from older trees were considerably greater than those from younger trees. Metabolomic analysis revealed that the concentrations of 18 amino acids and derivatives (including histidine and serine) in nuts from older trees were markedly higher than those in nuts from younger trees. Transcriptome and metabolomic correlation analysis identified 16 genes, including TgPK (pyruvate kinase), TgGAPDH (glyceraldehyde 3-phosphate dehydrogenase), and others, which exhibit higher expression levels in older trees compared to younger trees, as confirmed by qRT-PCR. These genes are associated with the biosynthesis of histidine, glutamic acid, tryptophan, and serine. Transient expression of TgPK in tobacco led to increased pyruvate kinase activity and amino acid content (histidine, tryptophan, and serine). Additionally, dual-luciferase assays and yeast one-hybrid results demonstrated that TgWRKY21 positively regulates TgPK expression by directly binding to the TgPK promoter. These findings not only demonstrate the nutritional differences between nuts from young and old trees but also offer fresh insights into the development of nutritional sources and functional components based on nuts from old trees, enriching our understanding of the potential benefits of utilizing nuts from older trees.

Peripatric speciation within Torreya fargesii (Taxaceae) in the Hengduan Mountains inferred from multi-loci phylogeography.

BACKGROUND: The Hengduan Mountains (HDM) are one of the major global biodiversity hotspots in the world. Several evolutionary scenarios, especially in-situ diversification, have been proposed to account for the high species richness of temperate plants. However, peripatric speciation, an important mode of allopatric speciation, has seldom been reported in this region. RESULTS: Here, two chloroplast DNA regions and 14 nuclear loci were sequenced for 112 individuals from 10 populations of Torreya fargesii var. fargesii and 63 individuals from 6 populations of T. fargesii var. yunnanensis. Population genetic analyses revealed that the two varieties are well differentiated genetically (FST, 0.5765) and have uneven genetic diversity (π, 0.00221 vs. 0.00073 on an average of nuclear loci). The gene genealogical relationship showed that T. fargesii var. yunnanensis is inferred as derived from T. fargesii var. fargesii, which was further supported by the coalescent simulations (DIYABC, fastsimcoal2 and IMa2). By the coalescent simulations, the divergence time (~ 2.50-3.65 Ma) and the weak gene flow between the two varieties were detected. The gene flow was asymmetrical and only occurred in later stages of divergence, which is caused by second contact due to the population expansion (~ 0.61 Ma) in T. fargesii var. fargesii. In addition, niche modeling indicated that the two varieties are differentiated geographically and ecologically and have unbalanced distribution range. CONCLUSIONS: Overall, T. fargesii var. fargesii is always parapatric with respect to T. fargesii var. yunnanensis, and the latter derived from the former in peripatry of the HDM following a colonization from central China during the late Pliocene. Our findings demonstrate that peripatric speciation following dispersal events may be an important evolutionary scenario for the formation of biodiversity hotspot of the HDM.

Torreya grandis oil attenuates cognitive impairment in scopolamine-induced mice.

The oil of Torreya grandis (TGO), a common nut in China, is considered to be a bioactive edible oil and has a great value in functional food development. In this study, the neuroprotective effects of TGO were investigated on a scopolamine (SCOP)-induced C57BL/6J mouse model. The mice were pretreated with TGO for 30 days (1000 mg per kg per day and 3000 mg per kg per day, i.g.). Behavioral tests showed that the supplementation of TGO could prevent the cognitive deficits induced by SCOP. TGO rebalanced the disorder of the cholinergic system by upgrading the level of acetylcholine. TGO also alleviated the over-activation of microglia and inhibited neuroinflammation and oxidative stress. Additionally, TGO could regulate the composition of gut microbiota, increase the production of short-chain fatty acids, and decrease the content of lipopolysaccharides in the serum. In conclusion, TGO has the potential to prevent loss of memory and impairment of cognition, which may be related to its regulation of the gut microbiota-metabolite-brain axis.

TgLCYB1 regulated by TgWRKY22 enhances the tolerance of Torreya grandis to waterlogging stress.

ß-Carotene functions in plant growth and development and plays an important role in resisting abiotic stress, such as drought and salt stress. The specific function and mechanism by which ß-carotene responds to waterlogging stress, however, remain elusive. In this study, we found that ß-carotene content and lycopene cyclase (TgLCYB1) expression, both in leaves and roots of Torreya grandis, were increased under waterlogging treatment. Subcellular localization assays indicated that TgLCYB1 was localized in the chloroplasts. Phenotypic, physiological, and metabolome analysis showed that overexpression of TgLCYB1 enhanced the tolerance of tomato plants to waterlogging stress. Furthermore, application of a LCYB enzyme inhibitor, 2-(4-chlorophenylthio)-triethylamine hydrochloride, markedly enhanced the sensitivity of T. grandis to waterlogging stress. In addition, yeast one-hybrid assay, the dual luciferase assay system, and real-time quantitative PCR indicated that waterlogging stress induced TgWRKY22 to increase TgLCYB1 expression by binding to the TgLCYB1 promoter. Collectively, our results indicated that TgWRKY22 positively regulated TgLCYB1 expression to improve the activities of antioxidant enzyme and increase the levels of some key metabolites, thereby relieving waterlogging-induced oxidative damage, and consequently modulating the waterlogging stress response. This study contributes to a more comprehensive understanding of carotenoid functions and the role LCYB genes play in plant stress response.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants.

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Transcription factors TgbHLH95 and TgbZIP44 cotarget terpene biosynthesis gene TgGPPS in Torreya grandis nuts.

Terpenes are volatile compounds responsible for aroma and the postharvest quality of commercially important xiangfei (Torreya grandis) nuts, and there is interest in understanding the regulation of their biosynthesis. Here, a transcriptomics analysis of xiangfei nuts after harvest identified 156 genes associated with the terpenoid metabolic pathway. A geranyl diphosphate (GPP) synthase (TgGPPS) involved in production of the monoterpene precursor GPP was targeted for functional characterization, and its transcript levels positively correlated with terpene levels. Furthermore, transient overexpression of TgGPPS in tobacco (Nicotiana tabacum) leaves or tomato (Solanum lycopersicum) fruit led to monoterpene accumulation. Analysis of differentially expressed transcription factors identified one basic helix-loop-helix protein (TgbHLH95) and one basic leucine zipper protein (TgbZIP44) as potential TgGPPS regulators. TgbHLH95 showed significant transactivation of the TgGPPS promoter, and its transient overexpression in tobacco leaves led to monoterpene accumulation, whereas TgbZIP44 directly bound to an ACGT-containing element in the TgGPPS promoter, as determined by yeast 1-hybrid test and electrophoretic mobility shift assay. Bimolecular fluorescence complementation, firefly luciferase complementation imaging, co-immunoprecipitation, and GST pull-down assays confirmed a direct protein-protein interaction between TgbHLH95 and TgbZIP44 in vivo and in vitro, and in combination these proteins induced the TgGPPS promoter up to 4.7-fold in transactivation assays. These results indicate that a TgbHLH95/TgbZIP44 complex activates the TgGPPS promoter and upregulates terpene biosynthesis in xiangfei nuts after harvest, thereby contributing to its aroma.

Ti3C2Tx MXene nanosheets enhance the tolerance of Torreya grandis to Pb stress.

As a 2D nanomaterial, MXene (Ti3C2Tx) has shown enormous potential for use in fields such as biomedical and environmental pollution. However, the utilization of MXene materials in plants has received little attention thus far. The efficient use of MXene materials in agriculture and forestry is first highlighted in this study. Phenotypic and physiological analyses indicated that MXene application significantly enhanced the tolerance of Torreya grandis to Pb stress by reducing Pb accumulation. Furthermore, we illustrated two independent mechanisms of MXene material in reducing Pb accumulation in T. grandis: 1) MXene converted the available form of Pb into stable forms via its strong Pb adsorption ability, resulting in a decrease of the available form of Pb in soils, and 2) MXene application obviously increased the cell wall pectin content to restrict more Pb in the cell wall by regulating the expression of pectin synthesis/metabolism-related genes (TgPLL2, TgPLL11, TgPG5, TgPG30, TgGAUT3 and TgGAUT12) in T. grandis roots. Overall, this finding provides insight into the application of MXene material in modern agriculture and forestry, which will facilitate the rapid development of nanotechnology in sustainable agriculture and forestry.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis.

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C18 Δ9-elongase and a C20 Δ5-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ5-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

A comprehensive metabolomics analysis of Torreya grandis nuts with the effective de-astringent treatment during the postharvest ripening stage.

Astringency removal is important for the quality of Torreya grandis nut and occurs after harvest. Here, we evaluated the effect of NaHCO3 treatment on astringency removal and compared the differential metabolites of the seed coat and kernel using a UHPLC QQQ-MS-based metabolomics approach. The result revealed the nut astringency was primarily enriched in the seed coat with more soluble tannins. The NaHCO3 treatment greatly shortened the de-astringency process, as indicated by a faster conversion of soluble tannins to insoluble tannins and more acetaldehyde production. Besides, a total of 293 metabolites, including 92 phenolic acids and 37 flavonoids, were tentatively characterized in the seed coat. A further comparative analysis of the metabolomics indicated epigallocatechin, gallocatechin, catechin, procyanidin B1, B2, B3 and C1 to be the major metabolites influenced by the NaHCO3 treatment. This study provides new insights regarding the metabolite differences of Torreya grandis nuts processed with different de-astringent treatments.

Ethylene treatment promotes umami taste-active amino acids accumulation of Torreya grandis nuts post-harvest by comparative chemical and transcript analyses.

Amino acids play critical roles in physiological processes and also contribute significantly to fruit quality. In this study, the effect of exogenous ethylene on amino acids metabolism and related genes expression in Torreya grandis were investigated. The results revealed that ethylene treatment (3000 µL L-1 for 24 h) significantly increased amino acids level. Umami amino acids were distinctly upregulated in ethylene-treated versus control nuts, with glutamic and aspartic acids to demonstrate 1.9-fold and 2.1-fold increase. Transcriptome analysis revealed that deferentially expressed genes were mainly enriched in alanine aspartate and glutamate metabolism. RT-qPCR confirmed that ethylene treatment up-regulated expression of their biosynthesis genes (TgGOGAT1, TgAATC1, TgAATC4) concurrent with suppression of their degradation enzymes (TgGS2, TgGAD1, TgGAD3, TgASNS1). Ethylene treatment appears to promote umami taste-active amino acids and improve T. grandis nut quality post-harvest.

Metabolomics reveal changes in flavor quality and bioactive components in post-ripening Torreya grandis nuts and the underlying mechanism.

Secondary metabolites are a group of small molecules with critical roles in plants fitness in addition to their potential bioactivities in humans. Most of these compounds are associated with the flavor and quality formation of fruits or nuts during the development or the postharvest stages. Change in metabolic profiles and shifts underpinning the post-ripening process in T. grandis nuts are not yet reported. In this study, a large scale untargeted metabolomics approach was employed in T. grandis nuts, revealing for a total of 140 differential accumulated metabolites. Among them, nearly 60% of metabolites belonging to terpenoids, coumarins and phenolic acids, and phytohormones were showed a gradual accumulation pattern, while most of compounds in flavonoids were decreased during post-ripening. An in-depth analysis of changes in these metabolite classes suggest a framework for post-ripening process effect associated with the postharvest quality of T. grandis nuts for the first time.

Identification of novel bioactive proteins and their produced oligopeptides from Torreya grandis nuts using proteomic based prediction.

Torreya grandis nut is a chief functional food in China consumed for centuries. Besides its rich protein composition, increasing studies are now focusing on T. grandis functional proteins that have not yet identified. In this study, liquid chromatography coupled with mass spectrometry detection of smaller and major proteins, revealed that the major peptide was 36935.00 Da. Proteome sequencing annotated 142 proteins in total. Bioactive proteins such as defensin 4 was annotated and its anti-microbial function was verified. Finally, functional oligopeptides were predicted by searching sequences of digested peptides in databases. Ten group of oligopeptides were suggested to exhibit antioxidant, Angiotensin-converting enzyme inhibition, anti-inflammatory. The predicted antioxidant activity was experimentally validated. It is interesting that a peptide GYCVSDNN digested from defensin 4 showed antioxidant activity. This study reports novel functional peptides from T. grandis nuts that have not been isolated and/or included as functional ingredients in nutraceuticals and in food industry.

In Vitro Antioxidant Properties and Phenolic Profile of Acid Aqueous Ethanol Extracts from Torreya grandis Seed Coat.

Torreya grandis is an important economic forestry product in China, whose seeds are often consumed as edible nuts, or used as raw materials for oil processing. To date, as an important by-product of Torreya grandis, comprehensive studies regarding the Torreya grandis seed coat phenolic composition are lacking, which greatly limits its in-depth use. Therefore, in the present study, the Torreya grandis seed coat was extracted by acid aqueous ethanol (TE), and NMR and UHPLC-MS were used to identify the major phenolics. Together with the already known phenolics including protocatechuic acid, catechin, epigallocatechin gallate, and epicatechin gallate, the unreported new compound 2-hydroxy-2-(4-hydroxyphenylethyl) malonic acid was discovered. The results of the antioxidant properties showed that both TE and 2-hydroxy-2-(4-hydroxyphenylethyl) malonic acid exhibited strong ABTS, DPPH, and hydroxyl radical-scavenging activity, and significantly improved the O/W emulsion's oxidation stability. These results indicate that the TE and 2-hydroxy-2-(4-hydroxyphenylethyl) malonic acid could possibly be used in the future to manufacture functional foods or bioactive ingredients. Moreover, further studies are also needed to evaluate the biological activity of TE and 2-hydroxy-2-(4-hydroxyphenylethyl) malonic acid to increase the added value of Torreya grandis by-products.

Spatial variation, sources identification and risk assessment of soil heavy metals in a typical Torreya grandis cv. Merrillii plantation region of southeastern China.

Torreya grandis (Torreya grandis cv. Merrillii) is a unique nut tree species in China. Currently, researches on Torreya grandis focus on nuts quality and yield, while few works are related to the soil quality of Torreya grandis plantation. In this study, the typical Torreya grandis production areas of Zhuji, Shengzhou, Keqiao and Dongyang cities along the Kuaiji Mountain were selected. A total of 121 topsoil samples (0-20 cm) were collected based on a grid of 1 km × 1 km. The results indicated that the average concentrations of Cd, Cr, Cu, As, Ni and Pb in soils were 0.12, 49.01, 27.95, 14.28, 26.97 and 40.28 mg kg-1, respectively. The concentrations of six heavy metals all exceeded the background values, and there were different degrees of pollution levels. The results of Moran's I indicated that the spatial high-high clusters of soil heavy metals were mainly distributed in Zhuji and the junction of Shengzhou and Keqiao. The partial least squares path analysis of structural equation modeling (PLS-SEM) showed that OM and soil nutrients had extremely significant effects on soil heavy metals. Sources identification of principle component analysis (PCA) and positive matrix factorization model (PMF) revealed that agricultural activities, natural factors and mining were the main sources of soil heavy metals. The human health risks caused by soil heavy metals pollution were generally acceptable based on Monte Carlo simulation method. For the heavy-metal polluted area, management measures should be considered in order to protect human health.